A  TECHNIQUE  FOR  IMAGE  MAGNIFICATION  USING 
PARTITIONED  ITERATIVE  FUNCTION  SYSTEM 


Suman  K.  Mitra,  C.A.  Murthy  and  Malay  K.  Kundu 

Technical  Report  99-08 


July  1999 


Center  for  Multivariate  Analysis 
417  Thomas  Building 
Penn  State  University 
University  Park,  PA  16802 


Research  work  of  authors  was  partially  supported  by  the  Army  Research  Office 
under  Grant  DAAH04-96-1-0082.  The  United  States  Government  is  authorized  to 
reproduce  and  distribute  reprints  for  governmental  purposes  notwithstanding  any 
copyright  notation  hereon. 


DTIC  QUAUry  maPBCJTBD  4 


93  MASTER  COPY 


KEEP  THIS  COPY  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approv9d 
0MB  NO.  0704^188 


rtoonma  oufottn  lor  tmt  eoiisaion  oi  lotowntiion  >t  to  avoroqo  i  noor  ot  roioonto.  #iaoo^  mo  ••no  lOf  mtifuaiono  toofcnm^  ^  -m. 

:nnor.oq  ano  mo  oolo  nooooq.  orvi  comoiotmq  arm  rovto^rmt)  mo  eoiioaion  oi  mformonon  Sono  commom  rogoram9  mia  Dwfooo  owonotoa  Qt  ay  QWy  MM  Oi  m« 

:-Moaion  ot  ^fofmaiion  .neiuomq  tu99oaiiona  tor  roouonq  mit  OufOon.  to  Waanmotoo  Hoooouanoro  SorvH^t  Qiroaoraio  w  mtorrna w  (^oiiona  w  Hopw.  J 


2.  REPORT  DATE 

3.  REPORT  TYPE  AND  OATES  COVERED 

July  1999 

Technical  -  99-08 

1  All 

1  -GcNCY  JSE  ONLY  (L0av0  Oiant^f 


4  -.TLE  AND  SUBTITLE 

A  Technique  for  Image  Magnification  Using  Partitioned 
Iterative  Function  System 


6  AoTHORiS) 

Suman  K.  Mltra,  C.A.  Mur thy  and  Malay  K.  Kundu 


DAAH04-96-1-0082 


7  ===IFOHMING  ORGANIZATION  NAMESIS)  AND  AODRESS(ES) 
Center  for  Multivariate  Analysis 
417  Thomas  Building 
Department  of  Statistics 
Penn  State  University 
University  Park.  PA  168Q2 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

99-08 


9  ^  =ONSGRING  -  MONITORING  AGENCY  M AMEi  S :  ANO  ADORESSlES) 

i  S.  Armv  Research  Office 

P  O.  Box  12211  _ _ 

Research  Triangle  Park,  NC  27709-221 1 


10  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


M 


o  3 


11  SUPPLEMENTARY  NOTES 

The  Views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  consinicd  as 
an  official  Depanment  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


I2a.  01STRI8UTION  /  AVAILABILITY  STATEMENT 


12  b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


1 3  abstract  ,  Maximum  200  worasi 

A  new  technique  for  image  magnification  using  the  theory  of  fractals  is  proposed.  The 
technique  is  designed  assuming  self  transformability  property  of  images.  In  particular, 
the  magnification  task  is  performed  using  the  fnictal  code  of  the  image  insteiul  of  the 
original  one  resulting  in  a  reduction  in  memory  requirement.  To  generate  the  fractal 
codes,  Genetic  Algoritlim  with  elitist  model  is  used  wiiich  greatly  decreases  the  .search 
for  finding  self  similarities  in  the  given  image.  The  article  presents  both  theory  and 
implementation  of  the  proposed  method.  A  simple  distortion  measure  .scheme  and  a 
similarity  criterion  based  on  just  noticeabhs  diircreuce  have  also  Ixieti  proposed  to  judge 
the  image  quality  of  the  magnified  image.  Comparison  with  one  of  the  mast  popular 
magnification  techniques,  the  nearest  neighbor  technique,  is  made. 


14  SUBJECT  TERMS  Image  Magnification,  Iterated  Function  System  (IFS) 

15.  NUMBER  IF  PAGES 

32 

ramexonea  i.Lci.cicj.vc 

16.  PRICE  CODE 

17  SECURITY  CLASSIFICATION 

:r  report 
U.NCLASSIFIED 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 

20.  limitation  OF  ABSTRACT 

UL 

NSN  7S40^)1  •280-5500 

StMioard  Fonn  aw  iHfv. 
pJSmM  ev  ANii  SM.  »18 

Plfio  36~S/S.5H-I^h 

A  Technique  For  Image  Magnification  Using 
Partitioned  Iterative  Function  System 

Suman  K.  Mitra,  C.  A.  Murthy  and  Malay  K.  Kundu 

Machine  Intelligence  Unit, 

Indian  Statistical  Institute, 

203,  B.  T.  Road,Calcutta  700035,  INDIA. 
Email:{res9432,murthy,malay}<§isical.ac.in 

Abstract 

A  new  technique  for  image  magnification  using  the  theory  of  fractals  is  proposed.  The 
technique  is  designed  assuming  self  transformability  property  of  images.  In  particular, 
the  magnification  task  is  performed  using  the  fractal  code  of  the  image  instead  of  the 
original  one  resulting  in  a  reduction  in  memory  requirement.  To  generate  the  fractal 
codes.  Genetic  Algorithm  with  elitist  model  is  used  which  greatly  decreases  the  .search 
for  finding  self  similarities  in  the  given  image.  The  article  presents  both  theory  and 
implementation  of  the  proposed  method.  A  simple  distortion  measure  .scheme  and  a, 
similarity  criterion  based  on  just  noticeable  difference  have  also  been  proposed  to  judge 
the  image  quality  of  the  magnified  image.  Comparison  with  one  of  the  most  popular 
magnification  techniques,  the  nearest  neighbor  technique,  is  made. 

Keywords:  Image  Magnification,  Iterated  Function  System  (IFS),  Partitioned  Itera¬ 
tive  Function  System  (PIFS),  Genetic  Algorithms  (GAs). 

1  Introduction  ' 

Image  magnification  ideally  is  a  process  which  virtually  increases  image  resolution 
in  order  to  highlight  implicit  information  present  in  the  original  image,  not  (wident 
as  such.  It  can  be  looked  upon  as  a  scale  transformation.  Image  magnification  is 
used  for  various  application  like  matching  of  images  captured  irsing  different  sfui.sors 
(having  different  capturing  resolutions),  satellite  image  analysis  [1,  2],  rnedicaJ  image 
display  etc.  Normally  the  image  is  represented  in  the  form  of  a  two  dimensional 
array  of  pixel  values  (matrix  form),  and  it  requires  large  memory  space.  The  memory 
requirement  for  storage  or  bandwidth  requirement  for  transmission  is  greatly  reduced 
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when  different  coding  schemes  are  used.  The  actual  requirement  (memory/bandwidth) 
is  dependent  on  the  size  of  the  image  and  the  method  used  for  coding.  Conventional 
magnification  operation  is  generally  performed  on  an  image  represented  in  the  form  of 
a  matrix  (normal  form).  Before  applying  magnification  technique,  any  coded  image 
has  to  be  converted  into  normal  form  through  decoding  process  which  requires  some 
computational  cost.  So  it  should  beneficial  if  the  magnification  could  be  done  during 
decoding  itself.  Moreover  bandwidth  requirement  of  an  image  transmission  system 
would  be  reduced  further  if  image  of  smaller  size  is  transmitted  but  at  the  r(!ceivirig 
end  a  magnified  version  is  generated.  With  this  problem  in  mind  an  attempt,  is  mad(^ 
to  propose  a  new  magnification  technique  which  can  be  applied  directly  on  the  cock^l 
version  of  the  image. 

Fractal  image  coding  technique  is  one  of  the  efficient  approximate  image  coding  tech¬ 
niques  currently  available.  In  image  coding,  the  reconstructed  image  produced  is  usu¬ 
ally  subjectively  very  close  to  the  original  image.  Actually  the  codes  of  an  image  must 
implicitly  carry  all  the  spatial  information  associated  with  the  image.  Besides  the  spa¬ 
tial  information,  the  fractal  codes  carry  the  information  of  the  self  similarities  present 
in  the  image.  These  self  similarity  property  is  also  exploited  in  the  proposed  image 
magnification  technique.  We  call  it  as  fractal  image  magnification  technique. 

Fractal  geometry  has  recently  come  into  the  limelight  due  to  its  uses  in  various  sci¬ 
entific  and  technological  applications,  specially  in  the  field  of  computer  based  image 
processing.  It  is  being  successfully  used  for  image  data  representation  [3,  4]  and  as 
image  processing  tool  [5,  6].  In  this  connection,  the  use  of  Iterative  Function  System 
(IFS)  and  Collage  theorem  [3]  have  shown  a  remarkable  improvement  in  the  quality  of 
processing  compared  to  that  obtained  using  existing  image  processing  techniques. 

A  fully  automated  fractal  image  compression  scheme,  known  as  Partitioned  Iterative 
Function  System  (PIFS)  of  digital  images  was  first  proposed  by  .lacquin  [7].  TIk; 
basic  idea  of  fractal  image  compression  or  to  find  the  fractal  codes  of  an  image  is  to 
approximate  small  blocks,  called  range  blocks,  of  the  given  image  from  large  blocks, 
called  domain  blocks,  of  the  same  image.  Thus,  to  find  the  fractal  codes  for  a  given 
image,  a  mathematical  transformation  for  each  range  block  is  to  found,  which,  when 
appropriate  domain  block  gives  rise  to  an  approximation  of  the  range  block.  This 
set  of  transformations,  obtained  by  partioning  the  whole  image  is  called  Partitioned 
Iterative  Function  System  (PIFS).  In  this  scheme  the  self  similarity  of  the  image  blocks 
are  obtained  locally  so  the  scheme  is  also  known  as  local  iterative  function  system  [8]. 

Several  researchers  have  suggested  different  algorithms  with  different  motivations  to 
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obtain  PIFS  for  a  given  image.  We  have  already  suggested  a  faster  algorithm,  to 
obtain  PIFS,  using  Genetic  Algorithms  (GAs)  [9, 10].  GAs  [11, 12, 13]  are  optimization 
algorithms  which  are  modeled  according  to  the  biological  evolutionary  processes.  The.se 
optimization  techniques  reduce  the  search  space  and  time  significantly. 

In  the  present  work  an  attempt  is  made  to  use  fractal  codes  as  an  input  to  a,ii  image 
magnification  system.  Some  of  the  popular  techniques  of  digital  magnification  of  im¬ 
ages  are  nearest  neighbor,  bilinear  and  bicubic  interpolations.  All  these  techniques  are 
based  on  surface  interpolation,  usually,  in  the  interpolation  techniques,  the  global  in¬ 
formation  is  often  ignored.  The  local  or  semi  global  information  is  generally  exploited. 
In  the  proposed  scheme,  the  magnification  task  has  been  performed  by  using  fractal 
codes  where  both  the  local  and  semi  global  information  are  used.  The  scheme  is  noth¬ 
ing  but  a  decoding  scheme  of  fractal  codes  which  gives  rise  to  a  magnified  version  of 
the  original  image.  The  article  reports  the  initial  results  of  magnification  by  a  factor 
which  is  multiple  of  two.  The  technique  uses  fractal  codes  which  are  obtained  by  a  GA 
based  technique  [9,  10].  Comparison  with  the  nearest  neighbor  image  magnification 
method  has  also  been  reported  here. 

In  the  magnification  techniques  the  distortion  due  to  blocking  which  is  a  local  phe¬ 
nomenon  is  very  usual.  To  quantify  the  amount  of  distortion,  the  widely  u.sed  distortion 
measure  is  the  mean  squared  error  (MSE)  or  some  other  form  of  it.  MSE  is  a  global 
measure  which  fails  to  account  properly  the  local  distortion  due  to  blocking.  But  the 
blocking  effects  are  very  much  sensitive  to  the  human  visual  system.  So,  to  quantify 
the  global  and  the  local  distortions  simultaneously,  a  new  distortion  measure  (fidelity 
criterion)  is  introduced. 

In  the  process  of  magnification,  the  magnified  image  should  be  visually  similar  to  the 
original  one.  Beside  the  visual  judgment,  we  have  proposed  here  a  similarity  criterion 
based  on  just  noticeable  difference  (JND).  As  the  sizes  of  the  magnified  image  and  the 
original  image  are  different,  the  similarity  between  them  can’t  be  judged  by  inspecting 
the  pixel  values  alone.  Hence  the  JND  based  scheme  is  proposed  in  this  regard. 

Theory  and  key  features  of  IFS,  magnification  using  PIFS  and  GA  are  outlined  in 
Section  2.  The  methodology  of  using  fractal  codes  for  magnification  of  a  given  image 
is  described  in  Section  3.  A  new  fidelity  criterion  to  judge  the  performance  of  the 
proposed  algorithm  is  discussed  in  Section  4.  JND  based  similarity  criterion  is  discussed 
in  Section  5.  Section  6  presents  implementation  and  the  results.  Discussion  and 
conclusions  are  provided  in  Section  7. 
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2  Theory  and  Basic  Principles 

The  detailed  mathematical  description  of  the  IFS  theory,  Collage  theorem  and  other 
relevant  results  are  available  in  [3,  14,  15].  Only  the  salient  features  are  discussed  here. 
The  theory  of  IFS  in  image  coding  and  PIFS  in  image  magnification  are  described  in 
the  following  subsections.  The  basic  principle  of  Genetic  Algorithms  is  also  described. 

2.1  Theoretical  Foundation  of  IFS 

Let  /  be  a  given  image  which  belongs  to  the  set  X.  Generally  X  is  taken  as  the 
collection  of  compact  sets.  Our  intention  is  to  find  a  set  T  of  affine  contractive  maps 
for  which  the  given  imago  /  is  an  approximate  fixed  point.  The  fixed  point  or  attractor 
“A”  of  the  set  of  maps  T  is  defined  as  follows  : 

limJF^(J)  =  A,  VJeA, 

and  -^(A)  =  A,  where  X^{J)  is  defined  as 

T^[J)  =  ),with 

x\J)  =  T{j),  yjex. 

Also  the  set  of  maps  T  is  defined  as  follows: 

d{J^[J\),T{J2))  <  s  d{Ji,J2)\  yj\,J2  ^  X  and  0  <  s  <  1.  (1) 

Here  “d”  is  called  the  distance  measure  and  “s”  is  called  the  contractivity  factor  of  T. 

Let  d(/,  <  c  (2) 

where  e  is  a  small  positive  quantity.  Now,  by  Collage  theorem  [1],  it  can  be  shown  that 

d(/,A)<^  (3) 

Here,  (A,  is  called  iterative  function  system  and  JF  is  called  the  set  of  fractal  codes 
for  the  given  image  /. 
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2.1.1  Image  Coding  Using  PIFS 


Let,  /  be  a  given  image  having  size  w  x  w  and  the  range  of  gray  level  values  be 
[0,<7].  Thus  the  given  image  /  is  a  subset  of  The  image  is  partitioned  into  n 
non  overlapping  squares  of  size,  say  6x6,  and  let  this  partition  be  represented  by 
^  ) ^n}-  Each  TZi  is  named  as  range  block  where,  n  =  j  x  f.  Let 

V  be  the  collection  of  all  possible  blocks  (within  the  image  support)  which  is  of  size 
26  X  26  and  let  X>  =  {T^i,X^2,  •  •  •  Each  Vj  is  named  as  domain  block  with 

m  =  {w  —  2b)  X  {w  —  26).  Let, 

=  {f  ■  i  /  is  an  affine  contractive  map}. 

Now,  for  a  given  range  block  7?.j,  let,  G  Tj  be  such  that 

rf(Ri,  h(Vi))  <  <m,  /(©,))  V/  €  Tj,  Vj. 


Now  let  k  be  such  that 

fi\k{'Dk))  =  mm{  (i(7^i,  fi\jiT>j))  }  (4) 

Also,  let  fi\k{T^k)  = 

Our  aim  is  to  find  7li\k  for  each  i  G  {1,2, •••,n}  or  in  other  words  for  each  range 
block  {7li)  we  are  to  find  appropriately  matched  domain  block  {Vk)  and  appropriately 
matched  map  (/t|fc).  Thus  Wj  =  {'Dk,fi\k}  is  called  fractal  code  for  TZi  and  the  set 
W  =  (Wi  ,  i  =  l(l)n}  is  called  the  PIFS  of  the  given  image  /.  Figure  1  illustrates 
the  mapping  of  domain  blocks  to  the  range  blocks. 


Figure  1:  Mapping  for  an  IFS  scheme 
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2.2  Image  Magnification  Using  PIFS 


The  affine  contractive  transformation  is  constructed  using  the  fact  that  the  gray 
values  of  the  range  block  arc  scaled  and  shifted  version  of  the  gray  values  of  doinaiu 
block.  The  transformation  fi\j,  defined  on  IR^  is  such  that  fi\j{Vj)  approximate  TZi. 
fi\j  consists  of  two  parts,  one  for  spatial  information  and  the  other  for  information 
of  gray  values.  The  first  part  indicates  which  pixel  of  the  range  block  corresponds  to 
which  pixel  of  domain  block.  The  second  part  is  to  find  the  scaling  and  shift  parameters 
for  the  set  of  pixel  values  of  the  domain  blocks  to  the  range  blocks. 

The  first  part  denotes  the  shuffling  the  pixel  positions  of  the  domain  block  and  can 
be  achieved  by  using  any  one  of  the  eight  transformations  (isometries)  on  the  domain 
blocks[7,  9].  Once  the  first  part  is  fixed,  second  part  is  estimation  of  a  set  of  values 
(gray  values)  of  range  blocks  from  the  set  of  values  of  the  transformed  domain  l)locks. 
These  estimates  can  be  obtained  by  using  the  least  square  analysis  of  two  sets  of  values 
[9,  10]. 

The  second  part  is  obtained  using  least  square  analysis  of  two  sets  of  gray  valiums  once 
the  first  part  is  fixed.  Moreover  the  size  of  the  domain  blocks  is  double  that  of  the  range 
blocks.  But,  the  least  square  (straight  line  fitting)  needs  point  to  point  corresiauKkuiee. 
To  overcome  this,  one  has  to  construct  contracted  domain  blocks  such  that  the  number 
of  pixels  in  the  contracted  domain  blocks  become  equal  to  that  of  range  blocks.  The 
contracted  domain  blocks  are  obtained  by  adopting  any  one  of  the  two  techniques.  In 
the  first  technique,  the  average  values  of  four  neighboring  pixel  values  in  a  domain 
blocks  are  considered  as  the  pixel  values  of  the  contracted  domain  blocks  [7].  In  the 
other  scheme,  contracted  domain  blocks  are  constructed  by  taking  pixel  values  from 
every  alternative  rows  and  columns  of  the  domain  blocks  [9,  10].  In  the  present  article 
we  have  adopted  the  first  one. 

Thus  fi^k  can  be  looked  upon  as  mixture  of  two  transformations,  fi\k  =  ti\kC  ,  whore, 
C  is  contraction  operation  and  tjn  is  transformation  for  rows,  columns  and  gray  values. 

n 

Here  we  have,  I  =  [jUi  and  using  (2)  we  have, 

■i=l 

n  ^  ^ 

d{  [j^i,  [jHiik  )  <  €  (5) 

i=  1  i=:  1 

Now,  let  Ad  be  a  magnification  operator  such  that 

d{  (6) 

1=1  1=1 
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where  ei  is  a  small  positive  quantity.  Now  by  (5)  and  (6)  we  have, 

n  n 

dHjTli,  [jMilluk)  )  <  e2  (7) 

2=1 

where  £2  is  a  small  positive  quantity.  Again,  we  have, 

~  —  ti\k^  i'^k)- 

So, 

n  n 

di  [j'Jli,  U  kk  C{Vk)  )  <  £2.  (8) 

1=1  i=l 

Now,  reconstruction  of  images  using  the  operator  M  should  be  an  inverse  of  contraction 
operation  using  the  operator  C.  So,  by  (7)  and  (8) 

d  (  M  {ni\k)  ,  ti\k{Vk)  )  <  £3.  (9) 

Hence,  by  (7), 

n  n 

d  (  >  U  kki'Dk) )  <  €4.  (10) 

i=l  i=\ 

Both  £3  and  £4  are  small  positive  quantities. 

Thus  from  (10),  it  is  clear  that  there  is  no  need  of  constructing  the  magnification 
operator  M,  only  the  second  part  of  the  fractal  codes  has  to  be  applied  on  the  domain 
block  to  get  an  image  which  is  very  close  to  the  given  image  /  and  this  image  has  size 
double  that  of  the  given  image. 


2.3  Genetic  Algorithms 

Genetic  Algorithms  (GAs)  [11,  12,  13]  are  highly  parallel  and  adaptive  search  and 
machine  learning  processes  based  on  a  natural  selection  mechanism  of  biological  genetic 
system.  Parallelism  of  GAs  depend  upon  the  machine  used  for  computations.  GAs  help 
to  find  the  global  optimal  solution  without  getting  stuck  at  local  optima  as  they  deal 
with  multiple  points  (called,  chrornosomes)  simultaneously.  To  solve  the  optimization 
problem,  GAs  start  with  the  chromosomal  (structural)  representation  of  a  parameter 
set.  The  parameter  set  is  coded  as  a  string  of  finite  length  called  chromosome  or  simply 
string.  Usually,  the  chromosomes  are  strings  of  O’s  and  I’s.  If  the  length  of  chromosome 
(string)  is  I  then  total  number  of  chromosomes  is  2h  To  find  a  near  optimal  solution, 
three  basic  genetic  operators,  i)Selection,  ii)Crossover  and  iii)Mutation  are  exploited 
in  GAs. 
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In  selection  procedure  the  objective  function  values  or  the  fitness  function  values  of 
each  individual  string  is  responsible  for  its  selection  as  a  new  string  in  the  next  mating 
pool.  We  have  used  the  elitist  model  of  GAs  where  the  worst  string  in  the  present 
generation  is  replaced  by  the  best  string  of  the  previous  generation.  The  crossover 
operation  on  a  pair  of  strings  emulates  the  mating  process  of  natural  genetic  .system. 
This  process  is  very  often  in  natural  genetic  system  and  thus,  a  high  probability  i.s 
assigned  to  indicate  the  occurrence  of  this  operation.  In  mutation  operation  every 
bit  of  every  string  is  replaced  by  the  reverse  character  (i.e.  0  by  1  and  1  by  0)  with 
some  probability.  Usually  a  low  probability  is  assigned  for  mutation  operation  and  tlu; 
occurrence  of  this  operation  is  guided  by  this  probability.  We  have  u.sed  tin;  varying 
mutation  probability  scheme  [16]  to  guide  the  mutation  operation  in  the  present  work. 
Starting  from  the  initial  population  (of  strings)  a  new  population  is  created  using 
three  genetic  operators  as  described  above.  This  entire  process  is  called  an  iteration. 
In  GAs  a  considerable  number  of  iterations  are  performed  to  find  the  optimal  solution. 
The  string  which  possesses  optimal  fitness  value  among  all  the  strings  is  called  the 
optimal  string.  The  optimality  of  the  fitness  value  of  strings  is  problem  dependent.  If 
the  problem  is  a  minimizing  problem,  the  lowest  fitness  value  is  taken  as  the  optimal 
one  and  the  maximum  fitness  value  is  selected  as  the  optimal  one  if  the  proljlem  i.s  a 
maximization  problem.  The  convergence  of  GAs  to  an  optimal  solution  is  assurc-.d  as 
the  number  of  iterations  increases  [17]. 

The  methodologies  to  obtain  magnified  images  from  fractal  codes  are  now  discussed 
below. 


3  Methodology 

So  far  we  have  discussed  how  to  apply  the  fractal  codes  or  IPS  to  get  a  magnified 
image  which  is  double  in  size  that  of  the  given  image.  On  successive  applications  of 
this  proposed  algorithm,  magnification  by  factor  4,  8,  16  etc.  can  also  be  achieved. 
But  the  first  task  is  to  obtain  the  fractal  codes  or  IPS  for  a  given  image. 

3.1  Construction  of  PIFS  for  magnification 

The  size  of  the  range  blocks  plays  an  important  role  in  image  compression  as  well  as 
magnification.  If  small  blocks  are  taken,  the  finer  details  of  the  image  are  preserved  and 
restored  in  the  decompressed  image  but  the  compression  ratio  will  be  less.  On  the  other 
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hand  more  compression  will  be  achieved,  at  the  cost  of  finer  image  details,  if  larger 
range  blocks  are  considered.  Thus  a  trade  olT  has  to  be  made  to  get  good  quality 
decompressed  image  as  well  as  considerable  amount  of  compression.  But  the  main 
task  in  magnification  is  only  to  restore  all  the  image  details  and  almost  no  emphasis  is 
given  on  the  amount  of  compression  achieved.  So,  in  this  case,  small  range  blocks  are 
considered  to  keep  track  with  every  minute  details  of  the  original  image. 

In  the  proposed  algorithm,  to  obtain  the  fractal  codes  of  small  range  blocks  of  a 
given  image,  the  blocks  are  first  classified  into  two  groups  using  a  simple  classification 
scheme  [10].  The  groups  are  formed  according  to  the  variability  of  the  pix(^l  values  in 
the  blocks.  If  the  variability  of  a  block  is  low  i.e.,  if  the  variance  of  tin;  ])ix(!l  valutas 
in  the  block  is  below  a  fixed  value,  called  threshold,  we  call  the  block  as  smooth  type 
range  block.  Otherwise  we  call  it  a  rough  type  range  block.  The  threshold  value  to 
separate  the  range  blocks  into  two  types  is  obtained  from  the  valley  in  the  histogram 
of  the  variances  of  pixel  values  of  the  blocks.  All  the  pixel  values  in  a  smooth  type 
range  block  are  replaced  by  the  mean  of  its  pixel  values.  So,  it  is  enough  to  store  only 
the  discretized  mean  value.  On  the  other  hand  for  each  rough  type  range  block,  the 
appropriately  matched  domain  block  as  well  as  appropriately  matched  transformation 
from  eight  possible  isometric  transformations  [7]  have  to  be  searched  out.  To  solve 
this  search  problem  a  GA  based  technique  [9,  10]  is  adopted.  GA  is  a  search  t(^chniquc 
which  finds  out  the  optimal  solution  faster  than  the  exhaustive  search  techniepie  [10]. 

3.2  GA  to  Find  PIFS 

The  parameters  which  are  to  be  searched  using  GA  are  location  (starting  row  and 
starting  column)  of  domain  block  and  its  eight  possible  isometric  transformations  [7]. 
The  realization  of  the  first  is  two  integer  values  between  1  and  w  —  2b  and  the  second 
can  take  any  value  between  1  and  8.  Binary  strings  of  length  I  are  introduced  to 
represent  the  parameter  set.  Here  I  is  chosen  in  such  a  way  that  the  set  of  2'-  binary 
strings  exhausts  the  .whole  parametric  space. 

A  string  indicates  the  location  and  the  isometric  transformation  of  a  domain  block.  In 
fractal  codes  we  are  to  find  an  appropriate  domain  block  and  an  appropriate  trans¬ 
formation  for  each  range  block.  In  other  words  we  are  to  find  the  appropriate  string 
for  each  range  block.  Out  of  2'  strings  a  few  strings  are  selected  randomly  to  start 
the  GA.  Starting  with  the  initial  mating  pool  and  using  the  three  basic  operations, 
new  populations  are  generated  in  each  iteration  of  the  GA.  After  a  large  number  of 
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iterations,  the  GA  will  produce  a  near  global  optimal  solution.  To  obtain  the  appro¬ 
priate  string  in  each  step  we  are  to  calculate  the  fitness  function  of  each  string  in  each 
iteration.  Mean  square  error  (MSE)  is  used  as  fitness  function  of  a  string.  In  each 
mating  pool,  the  strings  first  under  go  crossing  over  operation  pairwise  and  the  muta¬ 
tion  operation  is  applied  in  each  bit  of  each  string.  Using  the  fractal  codes  or  PIFS  by 
GA,  a  magnification  of  order  two  is  achieved. 

In  the  next  subsection  the  technique  for  successive  magnification  has  been  described. 


3.3  Successive  Magnification 


In  the  case  of  successive  applications  of  the  algorithm,  or  to  obtain  magnification  of 
order  more  than  two,  tlie  fractal  cochis  ikhhI  not  I)e  computed  afrcish.  Tlici  fractal  codes 
that  are  used  in  a  step  (of  order  which  is  multiple  of  two)  are  obtained  by  modifying 
the  fractal  codes  obtained  in  the  previous  step.  In  particular,  the  transformations  U|;t 
are  modified  by  using  the  image  that  is  already  obtained  in  the  previous  step.  The 
locations  of  the  appropriately  matched  domain  blocks  are  kept  fixed  in  all  the  steps. 
Only  the  size  of  the  domain  blocks  is  increased  in  the  modification  process.  Thus,  in 
particular,  only  the  gray  level  transformation  in  tjjfc  is  to  be  modified.  The  gray  level 
transformations  are  obtained  using  least  square  technique.  In  this  technique  a  straight 
line  is  fitted  with  two  sets  of  gray  level  values  of  which  one  is  from  range  block  and 
other  is  from  contracted  domain  block.  In  the  successive  magnification  scheme  these 
two  sets  are  enlarged.  These  enlarged  sets  are  divided  into  several  parts  and  separate 
straight  lines  are  fitted  for  each  parts  using  the  same  least  square  technique.  The  sets 
are  divided  into  4  parts  for  achieving  magnification  by  a  factor  4  and  divided  in  to  IG 
parts  for  the  magnification  by  a  factor  8,  and  so  on.  So,  the  number  of  fractal  codes, 
in  a  step,  becomes  4  times  larger  than  its  counter  part  in  the  previous  step.  So,  it 
is  enough  to  find  the  fractal  codes  in  the  case  of  magnification  by  a  factor  2  and  in 
other  cases  these  codes  are  modified  accordingly.  The  modified  codes  are  then  used  for 
achieving  magnification  greater  than  two. 

The  next  subsection  provides  a  discussion  on  the  magnification  by  any  order. 


3.4  Magnification  By  Any  Order 

In  this  article  fractal  image  magnification  algorithm  is  implemented  with  magnification 
factors  which  are  multiples  of  2.  But  in  practice  one  may  need  to  magnify  the  given 
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image  by  other  factors  too.  One  way  of  performing  the  magnification  of  order  k  is 
to  select  the  domain  blocks  which  are  k  times  larger  than  that  of  the  range  block. 
But  if  k  is  large  i.  e.  if  the  size  of  the  domain  blocks  are  very  large  compared  to  that 
of  range  blocks,  the  similarity  patterns  between  range  blocks  and  domain  blocks  will 
not  be  observed.  In  that  case  true  magnification  will  not  be  possible.  To  avoid  this, 
magnifications  of  order  4,  8  and  16  are  obtained  by  modifying  the  fractal  codes  of 
range  block  size  2x2  and  using  magnified  image  obtained  in  the  previous  step  {e.g. 
magnification  of  order  4  is  obtained  using  the  twice  magnified  image  and  the  fractal 
codes.)  Similarly  to  obtain  magnification  of  order  3  one  has  to  consider  the  domain 
blocks  which  are  three  times  larger  than  that  of  range  blocks.  Hence,  by  modifying 
these  codes,  magnification  of  order  which  are  multiple  of  three  can  be  obtaiiKsd. 

On  the  other  hand,  magnification  by  factors  which  arc  not  multiples  of  two  or  three 
can  also  be  achieved  by  eonsidcuing  the  normaliz(!d  distance  between  the  range  blocks 
and  their  matched  domain  blocks  in  the  original  image.  The  normalized  distance's  an; 
stored  in  the  PITS  codes  for  each  range  block.  Now,  to  achieve  magnification  of  ordew 
k,  from  PIFS  codes,  the  location  of  the  matched  domain  blocks  for  range  blocks  may 
be  obtained  by  multiplying  the  normalized  distances  by  k.  Once  the  matched  domain 
block  is  fixed  the  rest  is  same  as  magnification  by  factor  two.  While  performing  the 
magnification  task,  if  the  distance  between  the  range  block  and  the  domain  block 
appears  to  be  fractional,  one  has  to  discretize  it. 

The  PIFS  codes  provide  some  loss  of  information  in  the  reconstructed  image.  So,  the 
image  quality,  in  magnification,  will  be  decreased  with  the  increasing  value  of  k.  This 
problem  can  be  handled  by  splitting  the  magnification  of  k  factor  task  into  several 
steps  and  by  modifying  the  PIFS  codes  in  each  step.  Though  some  loss  would  be 
incurreded,  the  overall  gain  will  be  in  terms  of  storage  requirements  as  instead  of  the 
original  image  the  codes  are  used,  for  the  task  to  perform. 

In  the  next  section  we  shall  discuss  the  evaluation  criteria  to  judge  the  performance  of 
the  proposed  algorithm. 


4  Fidelity  Criterion 

It  is  necessary  to  judge  the  performance  of  the  proposed  fractal  based  image  magnifica¬ 
tion  algorithm  for  understanding  its  properties  and  fruitfulness.  For  this  purpose  one 
has  to  measure  the  distortion  between  the  given  image  and  the  reconstructed  image.  To 
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quantify  the  amount  of  distortion,  the  widely  used  distortion  measure  is  peak  signal  to 
noise  ratio  (PSNR)  which  is  a  function  of  Mean  Squared  Error  (MSE).  MSE  or  PSNR 
examines  the  similarities  between  two  images.  But  MSE  is  a  size  dependent  measure 
z.e.,  the  two  images,  under  consideration,  should  have  the  same  size.  Moreover  it  is  a 
global  measure  which  is  the  average  of  pixel  to  pixel  difference.  It  does  not  accurately 
indicate  the  large  and  significant  local  distortions  due  to  blocking  or  blurring  [20]  as  it 
deals  with  the  average  distortions.  But  the  blocking  effects  are  very  much  sensitive  to 
the  human  visual  system.  So,  one  has  to  think  of  a  size  independent  measure  which 
reflects  local  as  well  as  global  distortions  and  judges  the  performance  of  magnification 
methodology.  A  new  fidelity  criterion  whose  performance  is  also  similar  to  that  of  vi¬ 
sual  judgment  is  introduced  to  find  out  the  distortion  between  the  given  image  and  the 
magnified  reconstructed  image.  The  overall  performances  of  the  proposed  algorithm  is 
found  using  this  new  distortion  measure. 

4.1  Edge  Based  Distortion  Measure 


The  images  that  are  obtained  from  the  codes  usually  have  specific  artifacts  such  as 
blocking,  ringing  and  blurring.  Actually  these  artifacts  are  reflected  more  prominently 
in  the  high  frequency  component  of  the  image  and  are  very  sensible  to  tlu^  human 
visual  system  [18].  So,  in  our  proposed  distortion  measure,  we  have  tried  to  measure 
the  dissimilarities  in  edge  pattern.  For  simplicity  only  the  vertical  and  horizontal  edges 
have  been  considered. 

Both  the  images  are  first  partitioned  into  blocks  proportional  to  their  respective  sizes 
in  such  a  way  that  both  images  contain  equal  number  of  blocks.  The  error  is  then 
measured  block  wise  and  finally  the  average  error  is  computed.  To  detect  the  edges  of 
each  block  we  have  used  the  scheme  suggested  by  Ramamurthi  et  al.  [19].  The  edge 
blocks  consist  of  value  “0”  and  “1”  where,  “1”  represents  the  presence  of  edg(\  Now  it 
is  expected  that  the  original  and  the  magnified  blocks  should  have  same  type;  of  edg(! 
distributions.  In  other  words  the  expected  run  of  I’s  present  in  both  the  blocks  should 
be  same  if  normalized  by  their  respective  sizes.  Thus  the  distortion  measure  is  defined 
by  the  difference  between  the  normalized  “expected  run”  of  I’s  present  in  the  given 
image  and  in  the  magnified  reconstructed  image.  The  vertical  and  horizontal  edges  are 
considered  separately  and  then  averaged  to  give  rise  to  the  final  distortion  measure  of 
a  block.  The  algorithm  of  the  proposed  distortion  criterion  is  described  below. 
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4.1.1  Description  of  the  algorithm 


Step  1  :  Partition  the  images,  ly  and  I2  (with  size  of  h  less  than  size  of  h)  into 
square  blocks  such  that  the  number  of  partitions  is  same  in  both  the  images.  Let 
Pi  and  p2  (with  pi  <  P2)  be  the  sizes  of  the  square  blocks  for  the  images  /|  and  I2 
respectively.  Let  these  blocks  be  5ii,  B12,  •  •  ■ ,  and  j52i,  B22,  ■  •  • ,  B2n- 

Step  2  :  From  Bij,i  =  1,2  and  j  =  1,2,  •••,n  compute  gradient  matrices  or  the 
edge  images.  Let  G^j  and  be  respectively  the  horizontal  and  the  vertical  gradient 
matrix.  The  elements  of  the  gradient  matrices  are  all  either  0  or  1.  The  gradient 
matrices  are  defined  as  follows 

G'‘(m,n)  =  0  if  <  T 

2 

=  1  if  Otherwise. 

and 

Gij{m,n)  =  0  if  HZ:fr.+fj::+un  <  T 

2 

=  1  if  Otherwise. 

Here  =  Gray  level  value  of  (m,  n)th  pixel  in  a  block  and  T  =  A  prefixed  tlireshold 
value. 

Step  3  :  Find  the  expected  run  of  I’s  present  in  both  horizontal  and  vertical  dir(>,c;tion.s 
in  both  G^j  and  G'^j.  Let  L  be  the  random  variable  denoting  the  number  of  run  of  I’s 
in  a  particular  gradient  matrix  in  a  particular  direction.  Compute  , 

,  Elf{L).  Here  expected  run  of  I’s  is  defined  as 

Number  of  times  the  run  of  length  k  appears 
^  Total  number  of  runs  {of  all  possible  lengths)  present’ 


Now  compute 


^  +  Et’jL)  +  EgCL)  +  Etyq) 


Step  4  :  Normalize  Eij{L)  by  respective  block  size. 

u  »  =  l 

=  if  ,:  =  2 

Step  5  :  Compute  the  final  error  measure  E  between  the  two  given  images  B  and 
E  is  defined  as 

tWiAL)  -  Ey{L)}\ 

^  j  =  l 

The  next  section  discuses  the  JND  based  similarity  criterion. 
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5  JND  Based  Similarity  Criterion 


Just  noticeable  difference  (JND)  measures  the  amount  of  change  in  gray  value  of  a 
pixel  in  comparison  of  its  surrounding  pixels.  Usually,  JND  is  used  to  evaluate  the 
edges  present  in  an  image  [22].  Here  we  have  proposed  a  similarity  measure  based  on 
JND  to  judge  the  similarity  between  two  images  having  unequal  size.  In  particular, 
the  proposed  similarity  measure  judges  the  performance  of  the  proposed  PIFS  based 
image  magnification  technique.  JND  is  basically  the  difference  in  contrast  of  an  object 
from  its  background  and  it  plays  an  important  role  in  human  visual  system.  The 
human  visual  system  (HVS)  model  [20]  deals  mainly  with  three  factors,  the  luminance! 
level,  spatial  frequency  and  signal  content.  Out  of  these,  the  perceived  luminance  i.s 
a  nonlinear  function  of  incident  luminance.  According  to  Weber’s  Law  [2:1],  if  the 
luminance  {Lb  +  AL)  of  an  ol)ject  is  just  noticealdy  different  from  its  btickgrouiid 
luminance  Lb,  then  1\L/Lb  =  constant.  Therefore,  the  just  noticeable  difference 
(JND)  AT  increases  with  the  increase  in  L/j. 

In  the  present  case,  we  have  developed  a  criterion  based  on  JND,  to  judge  the  similarity 
between  two  images,  which  are  unequal  in  size.  In  particular,  first  the  average  values 
of  AL  from  both  the  images  are  computed.  Next  percentage  of  similarity  present  is 
computed  using  average  values  of  AL.  We  call  this  percentage  of  similarity  based  on 
JND  as  JND  similarity.  Note  that,  JND  similarity  100,  implies  complete  similarity 
between  two  images.  The  computation  of  Lb  and  AL  are  carried  out  as  described  in 
[22]. 

6  Implementation  and  Results 


To  find  the  fractal  codes  for  a  given  image  the  search  is  to  be  made  for  all  possible! 
domain  blocks  as  well  as  eight  possible  isometric  transformations  [7].  To  reduce  tlu! 
search  space  and  time  Genetic  Algorithm  is  used  instead  of  exhaustive  search.  The 
search  space  reduction  is  achieved  since  near  optimal  solutions  are  usually  satisfactory 
and,  intuitively,  solutions  which  are  far  away  from  the  expected  are  rejected  in  a 
probabilistic  manner.  This  is  the  reason  for  GA  to  perform  well  for  optimization 
problems.  The  good  performance  of  GAs  to  find  fractal  codes  of  a  given  image  has 
already  been  shown  by  Mitra  et.  al.[10].  The  results  are  quite  satisfactory  and  at  least 
20  times  reduction  in  the  search  space  is  achieved. 

For  the  specific  implementation  of  the  proposed  algorithm,  a  part  of  the  original  “Lena” 
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image  (Figure  2(a))  is  treated  as  the  original  input  image.  The  input  image  is  a 
128  X  128,  8  bit/pixel  image.  The  GA  based  technique  [10]  is  applied  to  generate 
the  fractal  codes.  Moreover,  a  simple  classification  scheme  [10]  for  range  blocks  have 
been  adopted  to  retain  the  regions  where  the  gray  level  variation  is  minimum.  In  the 
classification  scheme,  the  range  blocks  are  grouped  into  two  classes  viz.,  “smooth”  and 
“rough” .  Every  pixel  value  of  a  smooth  range  block  is  replaced  by  the  average  of  all 
the  pixel  values.  For  each  rough  type  range  block  the  GA  based  technique  [10]  is  used 
to  find  fractal  codes.  In  the  case  of  magnification  algorithm,  small  range  blocks  of  size 
2x2  are  considered  for  the  computation  of  fractal  codes.  It  is  true  that  compression 
ratio  will  be  reduced  by  considering  small  range  blocks  but,  the  finer  details  of  rough 
type  range  blocks  will  be  retained.  The  main  aim  of  a  magnification  task  is  to  magnify 
the  image  keeping  all  the  image  details.  So,  we  have  considered  small  range  blocks. 
Now,  Using  the  obtained  fractal  codes,  in  the  way  described  in  Section  3,  an  image  of 
size  256  x  256  is  reconstructed.  The  reconstructed  image  is  two  times  rnagnilied  than 
the  original  image.  This  image  is  found  to  be  very  close  to  the  original  image!  which  is 
judged  by  the  error  measure  and  the  similarity  measure  as  d(!,s(:ril)(!d  in  S(!c(.i()ii  A  and 
5  respectively.  The  fractal  codes  are  then  modified  stepwise,  as  described  in  Section 
3.3,  to  get  the  images  which  are  4  times  and  8  times  magnified  than  the  original  one. 
In  each  step,  the  error,  in  comparison  to  its  previous  step,  is  measured  successively. 
Also  the  similarities  of  magnified  images  are  judged,  successively,  by  the  JND  based 
similarity  criterion. 

The  proposed  algorithm  is  also  compared  with  the  nearest  neighbour  technique  for 
image  magnification  in  terms  of  proposed  distortion  measure  and  similarity  measure. 
Nearest  neighbour  is  the  simplest  method  of  digital  magnification.  Given  an  image  of 
size  w  X  w,  to  magnify  it  by  a  factor  k,  every  pixel  in  the  new  image  is  assigned  the 
gray  value  of  the  pixel  in  the  original  image  which  is  nearest  to  it.  This  is  equivalent 
to  repeating  the  gray  values  kx  k  times  to  obtain  the  magnified  image.  The  resultant 
image  for  large  magnification  factors  will  have  prominent  block  like  structures  due 
to  lack  of  smoothness.  The  other  techniques  of  digital  image  magnification  a,r(!  l^asi- 
cally  interpolation  methodologies  which  are  based  on  linear,  bilinear,  cubic  or  bicubic 
interpolation  [23,  24,  25,  26]. 

The  proposed  algorithm  has  also  been  implemented  on  a  “Low  Flying  Aircraft”  (LFA) 
image  having  size  128  x  128  and  range  of  gray  level  values  0  to  255.  Other  parameters 
of  the  algorithm  are  kept  fixed  as  in  the  case  of  “Lena”  image.  All  the  results  obtained 
are  presented  in  Table  1  and  Table  2. 


15 


Table  1:  The  results  obtained  in  terms  of  Distortion  of  the  Image  masnification  Algorithms 


Image 

— — - - - o  - 

Reconstructed  Image  Statistics 

MF 

Distortion 

MF 

Distortion 

MF 

Distortion 

Fractal 

NN 

Fractal 

NN 

Fractal 

NN 

Lena 

2 

1.18 

1.62 

4 

1.37 

3.11 

8 

1.24 

6.15 

LFA 

2 

2.32 

2.37 

4 

2.85 

4.59 

8 

2.43 

9.11 

MF=Magnification  Factor  and  NN=Nearest  Neighbour 


Table  2:  The  results  obtained  in  terms  of  similarity  of  the  Image  magnification  Algorithms 


Image 

Reconstructed  Image  Statistics 

MF 

Similarity  (%) 

MF 

Similarity  (%) 

MF 

Similarity  (%) 

Fractal 

NN 

Fractal 

NN 

Fractal 

NN 

Lena 

2 

63.47 

64.39 

4 

57.75 

49.06 

8 

54.50 

44.38 

LFA 

2 

82.05 

68.22 

4 

72.46 

51.83 

8 

64.40 

47.03 

MF=Magnification  Factor  and  NN=Nearest  Neighbour 


The  original  and  decoded  images  of  “Lena”  are  shown  in  figures  2  and  3  res]){H:tiv(;ly. 
Figure  4,  6  and  8  arc  respectiv(’.ly  the  two  times,  four  times  and  eight  tinuis  magnificHl 
images  of  “Lena”  using  the  proposed  fractal  based  technique  while  figures  5,  7  and 
9  are  respectively  two  times,  four  times  and  eight  times  magnified  images  of  “Lena” 
using  nearest  neighbour  technique.  Figures  10  and  11  are  respectively  the  original  and 
decoded  image  of  “LFA” .  The  results  of  fractal  based  magnification  of  “LFA”  image 
are  shown  in  figures  12,  13  and  14  while  figure  15  is  the  eight  times  magnified  image 
of  “LFA”  using  nearest  niegbour  technique. 

From  Table  1,  it  is  evident  that  in  terms  of  the  proposed  error  criterion  the  performance 
of  the  proposed  fractal  based  image  magnification  scheme  is  better  than  that  of  the 
nearest  nighbour  technique.  The  results  presented  in  Table  2  are  showing  that  in  terms 
of  the  similarity  criterion,  the  nearest  neighbour  technique  is  better  than  that  of  the 
fractal  based  technique  but  the  later  appears  to  be  better  for  magnification  factor  more 
than  two.  Comparing  figures  4  with  5,  6  with  7  and  8  with  9,  visually,  one  can  find 
some  ringing  and  blurring  are  present  in  the  case  of  nearest  nighbour  technique  for 
magnification  of  order  more  than  two.  On  the  other  hand,  in  the  case  of  proposed 
fractal  based  magnification  a  few  block  effects  have  been  observed. 
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Figure  2:  Original  “Lena”  image 


Figure  3:  Decoded  “Lena”  image 

7  Discussion  and  Conclusions 

The  most  important  advantage  of  the  proposed  technique  of  fractal  image  magnification 
is  that  it  utilizes  the  coded  (fractal)  version  of  the  input  image  instead  of  the  original 
image.  Therefore  it  is  cost  effective  in  the  sense  of  storage  space  and  time  as  no 
decoding  is  performed  at  the  receiving  end  in  case  of  transmission  of  the  codes. 

Another  advantage  of  fractal  image  magnification  is  that  it  magnifies  the  image  by 
expanding  the  fractal  codes  or  the  transformations  which  may  be  looked  upon  as 
independent  of  image  resolution.  The  only  error  involved  with  it  is  the  problem  of 
discretization.  Thus  the  structure  and  the  shape  of  the  image  remains  almost  same. 
In  a  sense,  it  is  like  interpolation  resulting  in  a  sharper  expanded  image.  Other  image 
magnification  schemes  use  pixel  replication  to  expand  image.  Pixel  replication  makes 
an  image  blocky,  blurry  and  patchy  after  a  certain  extent  of  expansion. 

The  size  of  the  range  blocks  considered  plays  a  vital  role  in  fractal  image  compression 
and  fractal  image  magnification  (ref.  section  3.1).  In  particular  these  two  algorithms 
are  in  opposite  direction  from  the  point  of  view  of  range  block  size.  So,  one  can 
think  of  an  optimal  range  block  size  for  which  good  quality  magnified  images  can  be 
reconstructed  from  the  fractal  codes  and  at  the  same  time  considerable  amount  of 
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Figure  5:  Two  times  magnified  “Lena”  using  Nearest  Neighbour 
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Figure  9:  Eight  times  magnified  “Lena”  using  Nearest  Neighbour 


Figure  10:  Original  “LFA”  image 


Figure  12:  Two  times  magnified  “LFA”  using  Fractal  technique 
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Figure  13:  Four  times  magnified  “LFA”  Fractal  technique 

compression  (in  terms  of  compression  ratio)  can  be  achieved.  To  solve  this  luoblem 
one  can  think  of  quadtree  partitioning  of  the  images  instead  of  square  partitioning 
while  generating  the  fractal  codes  [27].  Another  scheme,  to  obtain  the  PIFS  codes  of  an 
image,  has  been  suggested  by  Thomas  et  al.  [28]  can  also  be  adopted  in  this  connection. 
In  this  scheme  they  have  considered  irregular  shaped  range  blocks.  Automatically  the 
matched  domain  blocks  are  just  the  scaled  and  magnified  versions  of  these  irregular 
shaped  range  blocks. 

In  the  present  article  we  have  introduced  a  new  distortion  measure  or  fidelity  criterion 
to  judge  the  performance  of  the  proposed  algorithm.  There  are  other  methods  which 
are  non-parametric  statistical  tests  [29]  for  the  same  purpose.  The  common  tests 
for  examining  the  degree  of  association  between  two  distributions  whose  distribution 
functions  are  unknown  are  Sign  test,  Wald  Wolfowitz  Run  test,  Wilcoxson  test  and 
Kolmogorov  Smirnov  test.  Another  evaluating  criteria  based  on  fractal  dimension  has 
already  been  suggested  by  Lalitha  et  al.  [30].  But  the  most  important  feature  which 
should  be  considered  while  examining  the  distortion  between  two  images  is  tlie  edge 
distribution  of  the  images  as  the  edges  are  very  sensitive  to  human  eyes.  But  neither 
the  fractal  based  evaluating  criteria  nor  the  statistical  tests  take  care  of  distortions 
present  in  the  edges.  The  main  advantage  of  the  proposed  error  criterion  is  that  it 
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takes  care  of  distortions  in  the  edges.  Thus,  one  of  the  important  tasks  is  to  find  the 
proper  edges  in  the  images  for  the  implementation  of  the  proposed  distortion  measure. 
We  have  used  a  very  simple  technique  for  the  detection  of  edges  though  one  may 
suggest  more  complex  techniques  for  it. 

The  other  measure  proposed  for  judging  the  performance  of  the  proposed  fractal  based 
image  magnification  technique  is  JND  based  similarity  criterion.  This  measure  also 
takes  care  of  the  distribution  of  edges  as  JND  is  basically  the  change  in  contrast  of  an 
object  with  respect  to  its  background.  But  one  disadvantage  of  this  similarity  measure 
is  that  it  deals  with  the  change  in  pixel  values  ignoring  the  edge  pattern. 
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Summary 

Image  magnification  is  a  process  which  increases  image  size,  keeping  all  the  image 
details  unaffected,  in  order  to  highlight  implicit  information  present  in  the  image,  not 
evident  as  such.  Image  magnification  is  used  for  various  applications  like,  satellite 
image  analysis,  medical  image  display.  Generally,  the  irngc  is  represented  in  tlu;  form 
of  a  two  dimensional  array  of  pixels  values  (matrix  form)  and  it  requires  huge  memory 
space.  The  memory  requirement  of  storage  is  greatly  reduced  by  using  some  form  of 
image  compression.  But  for  performing  image  processing  tasks,  the  coded  image  is 
usually  brought  back  to  its  normal  form.  It  is  increasingly  the  case  that  the  coded 
form  of  the  image  is  used,  instead  of  the  normal  form,  as  the  input  to  perform  different 
image  processing  tasks.  With  this  aim  in  mind  an  attempt  is  made,  in  the  present 
article,  to  propose  a  new  magnification  technique  which  can  be  applied  directly  on  the 
coded  form  of  the  image.  In  particular  the  proposed  algorithm  is  using  that  coded 
form  of  the  image  which  is  obtained  by  fractal  image  compression  technique. 

Recently,  fractal  based  imago  c:()mprossioii  is  very  poi)ular  and  there  are  many  (,(h:Ii- 
niques  available  in  the  literature  for  finding  fractal  code  of  an  image.  The  encoding 
process  involves  in  computing  a  set  of  linear  contractive  maps  from  the  target  image. 
In  the  decoding  process,  the  obtained  set  of  maps  is  applied  on  an  arbitrary  image  in 
an  iterative  way  to  result  in  an  image  which  is  very  close  to  the  target  image.  The  set 
of  maps  is  called  fractal  code  or  Partitioned  Iterative  Function  System  (PIFS)  code  of 
the  image.  In  the  process  of  iterative  sequence,  PIFS  code  converges  to  a  fixed  image 
which  is  very  close  to  the  target  image.  Computational  task  for  finding  PIFS  code  of 
an  image  is  usually  time  consuming.  But  we  have  used  here  a  cost  effective  Genetic 
Algorithm  (GA)  based  technique  to  find  PIFS  code  to  propose  a  new  magnification 
technique. 

Conventional  magnification  techniques  are  basically  interpolation  methodologies  which 
are  based  on  linear,  bilinear,  cubic  or  bicubic  interpolations.  In  the  proposed  algorithm 
magnified  version  of  an  image  is  obtained  tising  the  reconstruction  of  the  frael.al  code; 
of  that  image.  No  magnification  operator  like  interpolation  is  needed.  Only  the  fractal 
code  or  in  otherwords  the  set  of  affine  contractive  maps  is  needed  to  magnify  the  image. 
We  have  proposed  here  the  technique  of  magnification  of  orders  which  are  multiple  of 
two.  The  technique  can  be  extended  to  the  case  of  magnification  by  any  order. 

We  have  also  proposed  two  techniques  to  judge  the  performance  of  the  proposed  magni¬ 
fication  algorithm.  The  main  task  involves  here  to  measure  the  distortion  or  similarity 
between  the  given  image  and  the  magnified  image.  The  commonly  used  distortion  rnea- 
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sure  is  peak  signal  to  noise  ratio  (PSNR)  which  is  a  function  of  mean  squared  error 
(MSB).  The  MSB  is  a  global  measure  and  also  it  is  image  size  dependent.  The  pro¬ 
posed  techniques  are  not  only  image  size  independent  but  also  utilizes  both  global  and 
local  information.  The  first  technique  is  a  distortion  measure  based  on  the  edge  distri¬ 
bution  of  the  images  and  indicates  the  influence  of  artifacts  like  blocking,  blurring  and 
ringing  which  may  appear  due  to  magnification.  The  other  one  is  a  similarity  measure 
based  on  just  noticeable  difference  (JND)  which  is  nothing  but  change  in  luminance  of 
an  object  pixel  with  respect  to  its  background  pixels.  The  overall  performance  of  the 
proposed  magnification  technique  is  found  to  be  satisfactory  both  qualitatively  and 
quantitatively.  Compjui.sou  with  one  of  tin;  inost  j)()])ul}u'  ina.giiifica.tioii  t(!c,liiii(|ii(;s, 
the  nearest  neighbor  technique,  is  made. 
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